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Radio and X-Ray Detectability of Buoyant Radio Plasma 
Bubbles in Clusters of Galaxies 
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Abstract. The Chandra X-ray Observatory is finding a surprisingly large number of cavities in the X-ray emitting 
intracluster medium (ICM) , produced by the release of radio plasma from active galactic nuclei. In this Letter, we 
present simple analytic formulae for the evolution of the X-ray deficit and for the radio spectrum of a buoyantly 
rising bubble. The aim of this work is to provide a theoretical framework for the planning and the analysis of X-ray 
and radio observations of galaxy clusters. We show that the cluster volume tested for the presence of cavities by 
X-ray observations is a strongly rising function of the sensitivity. 
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1. Introduction 



An early key result from the Chandra X-ray observatory 
was the discovery of numerous X-ray cavities in clusters of 
galaxies, preceded by pioneering detections by the ROSAT 
satellite. For example, cavities were found in the Perseus 
cluster (Bohringer et al. 1993), the Cygnus-A cluster 
(Carilh et al. 1994; Fabian et al. 2000), the Hydra-A clus- 
ter (McNamara et al. 2000), AbeU 2597 (McNamara et al. 

2001) , AbeU 4059 (Huang & Sarazin 1998; Heinz et al. 

2002) , AbeU 2199 (Fabian 2001), AbeU 2052 (Blanton 
et al. 2001), close to M84 in the Virgo Cluster (Finoguenov 
& Jones 2001), in the RBS797 cluster (Schindler et al. 
2001), and in the MKW3s cluster (Mazzotta et al. 2001). 
In most of these cases, the cavities are clearly coincident 
with the lobes of a radio galaxy at the cluster center. 
However, some clusters exhibit also cavities without de- 
tectable radio emission, namely in Perseus, AbeU 2597, 
and AbeU 4059. The latter class of cavities are also be- 
lieved to be filled with radio plasma, but during the buoy- 
ant rise of the very light radio plasma in the cluster at- 
mosphere (Gull & Northover 1973; Churazov et al. 2000, 
2001; Briiggen & Kaiser 2001; Churazov et al. 2002) adi- 
abatic expansion and synchrotron/inverse Compton radi- 
ation losses should have diminished the observable radio 
emitting electron population, leaving behind a so called 
ghost cavity or radio ghost (Enfilin 1999). The detectabil- 
ity of an X-ray cavity should also decrease with increasing 
distance from the cluster center due to the decreasing ratio 
of the missing X-ray emission from the volume occupied 
by the bubble to the fore- and background X-ray emission. 



2. Rising buoyant bubbles 

While the early phase of radio galaxy evolution is char- 
acterized by supersonic expansion into the surrounding 
medium, the radio lobes quickly settle into pressure equi- 
librium with the ICM after the AGN has shut off. Our 
description sets in at this moment ii, where the bubble 
is located at a cluster radius ri with volume Vb.i and 
pressure Pi = PicM{ri)- The bubble will then quickly 
approach a terminal velocity Vh{r)^ governed by the bal- 
ance of buoyancy and drag forces. During its rise, the 
bubble volume changes according to the adiabatic law 
Vh{r) — Vb,i (P(r)/Pi)^^/^'p, where the adiabatic index 
7rp is close to 4/3, which we will take as our fiducial 
value for numerical examples. The magnetic field strength 
should evolve according to B{t) = Bi (Vb(i)/Vb,i)^^^^ = 
Bi (P(t)/Pi)2/(3T'p), if the expansion of the bubble is 
isotropic. For simplicity, we only consider spherical bub- 
bles with radius rb, which gives sufficiently accurate esti- 
mates for most applications. If the bubble becomes highly 
deformed or even disintegrates, more sophisticated mod- 
els than ours will have to be used. We further assume, 
that entrainment of environmental gas into the bubble is 
dynamically insignificant on the considered time-scales, 
implying that the bubble is X-ray dark. Numerical simu- 
lations (e.g., Reynolds et al. 2001) support the latter. 

It is often convenient to express physical quantities 
like the bubble volume and the magnetic field strength 
in terms of the values they would have if the bubble 
were adiabatically moved to the cluster center. We de- 
note these by the subscript (e.g., Vh.o — Vh{r = 0) = 
Vb,! (Pn/Pi)^^^'^''')- As our working example, we will in- 
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vestigate a cluster described by an isothermal /3-profile 
with a density profile of p(r) = po (l + (''/?'c)^)~^^^^, pres- 
sure P{r) = Po (1 + {r/rc)'^)~^^/'^, and constant sound 
speed Cs(5. 

We define the origin of the cluster coordinate system at 
the cluster center, with the x and y axis defining the image 
plane and the x-axis the line of sight to the observer, and 
the coordinate system oriented so that the bubble center is 
located in the x-z plane at r = (r cos^, 0, r sinO). Its pro- 
jected distance from the cluster center is i? = fxr, where 
H — cos 9. The angle should be roughly conserved along 
the bubble's trajectory in a spherical cluster atmosphere. 

The buoyancy speed of the bubble can be estimated 
from the balance of buoyancy and drag forces. The buoy- 
ancy force 



dP 



buoyancy — Tj, 5 p — 7^''^ 



(1) 



depends on the gravitational acceleration g, which can be 
expressed via the pressure gradient in a cluster with hydro- 
statical equilibrium. The hydro dynamical drag of subsonic 
motion is well approximated by 



-Fdrag = CdTrrbPV^, 



(2) 



with Cd « 0.5. Using rb = rb,o(P/Po)"^/^^^'p\ the bubble 
velocity is given by 




4rb dP 



3 Cd p dr 
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(3) 



(4) 



is the characteristic buoyancy timescale across one core 
radius. The typical rise velocity is a fraction (oc ^Jr^Jrc) 
of the sound speed Cg. The time the bubble needs to travel 
from ri to r2 is then 



i(ri,r2) = 



B 



4' 47rp 



(5) 



where Ba;(a, b) = dyy"" {l — y f denotes the incomplete 
beta function, and [f{x)\%\ = f{x2) — f{xi). 

If the initial bubble is large compared to the cluster 
core Eq. 3 can give supersonic rise velocities. In such a case 
Eq. 2 is no longer valid. Instead, strongly increased dissi- 
pation will limit the rise velocity to the subsonic regime. 
In such a case one will adopt Wb = cthC-s, with Ob ~ 0.5. 

3. X-ray deficit evolution 

The X-ray emissivity at cluster radius r and for a given 
density and pressure profile is = K^{T{r)) n{r)'^ , where 
n is the electron number density and Ay(T) is the plasma 
cooling function. 

In the following we will assume that the gas surround- 
ing the bubble has essentially settled back into hydrostatic 



equilibrium in the dark matter potential, which is unper- 
turbed by the radio galaxy activity. 

In general, the surface brightness integral must be 
solved numerically, but in the case of a /3-modcl atmo- 
sphere, it can be expressed in closed form. For an undis- 
turbed line of sight (i.e., not intersecting a bubble) the 
observed surface brightness (neglecting cosmological cor- 
rections) is given by 



Ix,ciix,y) = lo ( 1 + 



y 



(6) 



where Iq = At ^c (47r)~^ B(3/3 — 5,5) is the central, 
and 7bg the background photon flux. B(a, b) is the beta 
function Bx=i{a, b). 

For a line of sight intersecting the surface of a bubble, 
the intersection points are located at {x,y, z±), with 



z±=r 1/1 - P^ ± 



(x — r iiY- 



(7) 



The surface brightness in the region occupied by the bub- 
ble is given by 

I^,h{x, y) = I^A^, y) - /o 1 + (8) 



sgn(z) B^(,) (1,3/3-1) 
2 B(3/3-ii) 



where we defined w{z) = z / {r^ + x +y + z ). 

Another important diagnostic is the number of missing 
counts ATVx from the cavities, i.e., the photons emitted 
from a spherical region of radius rb at distance r from the 
center in an unperturbed cluster atmosphere during the 
exposure interval fobs- For a detector with effective area 
Ai, and a source distance of D we define 



Then, 



dv 



2^3 



6/3-2 



o 1 ,,2-30 „ 



2)r, 



(9) 



(10) 



where v± = 1 + [(rb ± r)/rc]^. 

In order to design searches for cavities in promising 
cluster candidates, it is useful to estimate the expected sig- 
nal to noise ratio {S/J\f) of a potential observation. The 
signal is the expected number of missing photons from 
the volume occupied by the bubble, S = ANx- The back- 
ground, on top of which the missing photons have to be 
detected, is the expected number of photons A^exp from the 
area on the sky subtended by the cavity, assuming there 
were no cavity. This assumes that the instrument point 
spread function is not more extended than the cavity. The 
noise is then given by the fluctuation in this photon num- 
ber, which is A/' w 

The signal-to-noise ratio of a small bubble (so that 
the cluster density does not vary significantly across its 
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diameter, <C Max[rc,r]), in an isothermal /3-model is 
well approximated by 



4r2,ono(l+rVr2)-3''+-r 



:B(3/3- 



(11) 



where S is the ratio of the uniform background to central 
cluster photon flux. In most applications ^ <C 1. We write 



= ^(-0) = 



4 no 



b,0 



TT Xj 



3^/^T^ VrcB(3/3-i,i) 



(12) 



for the signal to noise ratio of a cavity moved to the clus- 
ter center. We write {S/Af)jj^iji for the minimum signal 
to noise required for a firm bubble detection. We define 
X = (5/7V)min/(5/A/')o, which is < 1 for all detectable 
bubbles. The subvolume of the cluster within which a bub- 
ble of the chosen size can be detected is shown in Fig. 
1. The dependence of the total integrated cluster volume 
Vdet, at which such a bubble can be detected, is displayed 
in Fig. 2 as a function of X and S. For 6 = the volume 
integral can be estimated analytically, and yields 



(13) 



where F = 1 - X«, ^ = 1/(3/3 - /3/7rp), = 6/(6/3 + 1 - 
4^/7,p), and r? = (6/3 - l)/(6/3 + 1 - 4/3/7,p). 

If the question of interest is the detectability of cavities 
of a given fixed size at all cluster radii one can use the 
incompressible limit 7rp ^ oo of Eqs. 11 and 13. This 
case is also displayed in Figs. 1 and 2. 

4. Radio spectrum evolution 

Synchrotron, inverse Compton, and adiabatic losses cool 
an electron with initial dimensionless momentum p\ = 
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Fig. 1. Contours of the X-ray deficit significance [S /N) of a 
bubble in a galaxy cluster. The contours mark locations at 
which a bubble has a significance which is lower by powers of 2 
than its significance if located at the cluster center (in the lower 
left corner of each figure, vertical axis is parallel to the line of 
sight). Left: the bubble volume expands adiabatically (with 
7rp = 4/3) with the pressure of the isothermal cluster ((3 = 
3/4). Right: the bubble volume is assumed to be independent 
of location (7rp oo). In both figures, an X-ray background 
with 1/100 of the central cluster surface brightness is assumed 
{6 = 0.01). 
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Fig. 2. Volume of the cluster [in units of r'c] within wliicli the 
signal to noise ratio {S/J\f) of the X-ray deficiency of a radio 
bubble exceeds a minimum, for the detection required value 
(<S/7V)min- The latter is given in units of {S/J^)i,), the signal 
to noise of a comparable bubble moved to the cluster center. 
The solid lines are for an adiabatically expanding bubble (7rp = 
4/3), and the dashed lines for an incompressible bubble (7rp — » 
oo). The background to central cluster X-ray surface brightness 
ratio takes the values 5 = 0, 0.001, 0.01, and 0.1 from top to 
bottom. The (top) curves with 6 = are given by Eq. 13. We 
adopted f3 = 3/4. 



momentum/ (mo c) down to p{pi,t) = Pi/{piq{t) + 
iVb{t)/Vb,i)i), where 



Pn 



ao 



n UB{r)+uc ( P{r) \ 



3 7rp 



(14) 



is the inverse maximal momentum. Here, oq = 
|(TT/(™ec). The magnetic and photon (mostly CMB) 
energy densities are denoted by Ub = B^/(87r) = 
UB,a (P/ Pq)^^^'^'^'''^'' and uc respectively. We have re- 
written the time integral in EnBlin & Gopal-Krishna 
(2001) as an integral over the trajectory of the bubble from 
ri to r2 . For our isothermal cluster this integral evaluates 
to 



wb,oB ^2 (k, Cb)+'McB ^2 (k, Cc) 



r2/r 



ri/r 



(15) 



with Cb = (11/3 - 37rp)/(47rp), CC = (3/3 - 37rp)/(47rp), 

K = 1/4, and r given by Eq. 4. 

If a bubble is rising with a constant velocity vy, (e.g. 
because it approaches the sound velocity) the following pa- 
rameters give the correct q value: Cb = (5/3 — 7rp)/(27rp), 
Cc = (/? - 7rp)/(27rp), K = 1/2, and t = Tc/ub- 

An initially relativistic power-law electron spectrum 
f{p, ti)dp = foP~^ dp for pmini <p< Pmaxi bccomes 



f{p,t) = fo{P{t)/Pi)^p-' {l-pq{t)) 



s-2 



(16) 
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Fig. 3. Bubble's central X-ray contrast (compared to the 
undisturbed cluster) for various angles between plane of sky 
and Bubble's trajectory, its radio flux, and its rising time as a 
function of the (unprojected) radial position. Adopted param- 
eters: rc = 20 kpc, Cs = 1400 km/s, /3 = 3/4, n = 10 kpc, 
rb,i = 8 kpc. Bo = lO/iG, 5 = 10-^ s = 2.4, Fi.4GHz(ri) = 
0.35 Jy. 

for Pmin{t) = p{Pminl,t) <P< Pmax(0 = P(Pmaxl, The 

synchrotron luminosity of these electrons is 

L^{t)= / dpP,{p,B)f{p,t), (17) 
Jo 

where Pv{p, B) is the synchrotron kernel. In the following 
we use the monochromatic approximation, which means 
that the total synchrotron radiation losses of an electron 
{Eg(p) = —qqUbP'^) are emitted at a single characteris- 
tic frequency v(p,B) — Ag B p^ with Ag = 3e/(27rTOoc). 
In this approximation, the synchrotron spectrum extends 
up to t-max = AsB2Pmax(*) < KB2q~'^{t) and can be 
written analytically as 

where 

_ 5-1 

was the original luminosity at r\. Eq. 18 (or the more 

accurate Eq. 17), and Eq. 15 (or Eq. 14) can be used to 
calculate the optical thin part of the radio spectrum of the 
source. 

5. Conclusion 

As an aid and stimulus for future observations we pro- 
vided simple analytic formulae of X-ray and radio prop- 
erties of rising buoyant bubbles of radio plasma in galaxy 
clusters. The detectability of a bubble decreases with its 
cluster radius. The X-ray contrast of a bubble moving in 



the plane of the sky decreases slowly until the X-ray back- 
ground dominates, then it drops rapidly. The contrast of a 
bubble moving along the line of sight declines quickly out- 
side the cluster core. Similarly, the radio luminosity of the 
bubble at a given frequency declines rapidly with increas- 
ing cluster radius and suddenly vanishes when cooling has 
removed the emitting electrons. After that point, only a 
weak flux of synchrotron-self Comptonized emission re- 
mains (EnBlin & Sunyaev 2001). Fig. 3 illustrates these 
dependencies for an example with parameters similar to 
Perseus A. 
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